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Abstract. For many classically chaotic systems, it is believed 
that in the semiclassical limit, the matrix elements of smooth ob- 
servables approach the phase space average of the observable. In 
the approach to the limit the matrix elements can fluctuate around 
this average. Here we study the variance of these fluctuations, for 
the quantum cat map on T 4 . We show that for certain maps and 
observables, the variance has a different rate of decay, than is ex- 
pected for generic chaotic systems. 



1. Introduction 

The basic question of "Quantum Chaos" is what signatures of the 
classical (chaotic) dynamics one can still find in the quantized system. 
One important such signature is given by the asymptotic behavior of 
expectation values, and the limit distribution around their average in 
the semiclassical limit. 

For classically chaotic systems, the average of smooth observables on 
phase space, along the trajectories of the system, tend almost always 
to the phase space average (i.e., the classical dynamics is ergodic). The 
quantum analog of this, following the correspondence principle, is that 
the expectation values of an observable should tend to the phase space 
average of the observable in the semiclassical limit. The first result 
in this direction is "Shnirel'man's theorem" , which states that at least 
in the sense of convergence in the mean square, when taking Planck's 
constant h — > 0, the matrix elements converge to the phase space av- 
erage j2 El UH 120] • This notion is usually referred to as "Quantum 
Ergodicity" (Q.E.). 1 

In the approach to the limit, the different matrix elements can fluc- 
tuate about their classical limit. In [7j Feingold and Peres proposed a 
formula for the variance of the fluctuations in the semiclassical limit. 
As mentioned above, according to Shnirel'man's theorem the variance 

1 The case where this convergence holds for all eigenstates is referred to as "Quan- 
tum Unique Ergodicity" (Q.U.E.). 
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vanishes as h — ■> 0. Nevertheless, after normalizing by an appropriate 
power of h, the Feingold-Peres formula relates the variance of the ma- 
trix elements to classical correlations of the observable. In Eckhardt 
et al. developed a semiclassical theory for the variance of these fluc- 
tuations, giving support for the validity of the Feingold-Peres formula, 
and suggesting that after normalizing by the variance, the fluctuation 
should be Gaussian (at least for hyperbolic systems). The analysis in 
[5], for the fluctuations of the matrix elements, can be extended to 
deal with quantum maps ^Bj- In particular, for generic hyperbolic 
quantum maps on the 2d-torus T 2d , it predicts Gaussian distribution 
with variance of order ™, where N — 1/h plays the role of the inverse 
of Planck's constant, and a is related to the classical variance of the 
observable. 

A fundamental example for a quantum map on the torus is the quan- 
tum cat map, originally introduced by Hannay and Berry 9j. Here the 
classical dynamics is simply the iteration of a symplectic linear map 
A e Sp(2d, Z) on T 2d = R 2d /Z 2d . When A has no roots of unity for 
eigenvalues, the classical dynamics is 'chaotic' (i.e., ergodic and mix- 
ing). 

For the quantization of the cat map, the admissible values of Planck's 
constant are inverses of integers, and the space of quantum states is 
then T~Cn = L 2 [(Zj/NZ*) d ]. The semiclassical limit is achieved by taking 
N — > oo, and we restrict the discussion to the case where N is prime. 
For / G C°°(T 2d ) a smooth observable, we denote by Op 7V (/) : Hn — ► 
its quantization. The quantization of the map A, is then a unitary 
operator U^{A) acting on Ti N . 

For d = 1, this model was studied extensively and it exhibits some 
interesting features El E3 EH- It is shown to be (Q.E.), but not 
(Q.U.E.). 2 However, in jT^j Kurlberg and Rudnick introduced a group 
of symmetries of the system, a family of commuting unitary maps of Ti^ 
that commute with Un(A). These operators are called Hecke operators, 
in analogy to a similar setup on the modular surface JUj- The space 
Hn has an orthonormal basis consisting of joint eigenstates {<^-}jLi 
called "Hecke eigenstates" . For these states the system is shown to be 
(Q.U.E.) 

In ^H] Kurlberg and Rudnick studied the fluctuation of the ma- 
trix elements for the desymmetrized system and gave a conjecture for 
the limit distribution, which is radically different from the behavior 



2 Faure, Nonnenmacher and De Bievre constructed a sparse subsequence of eigen- 
states where the expectation values do not converge to the phase space average ■ 
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expected in the generic case. In particular, a fourth moment calcula- 
tion for the distribution showed that it is not Gaussian. Furthermore, 
though the variance of the fluctuation is of order as expected, the 
classical factor is different from the classical variance, in contrast to 
the Feingold- Peres formula (for comparison see [To]). 



Remark 1.1. In [17 j Luo and Sarnak showed similar behavior of the 

variance for Hecke eigenforms on the modular surface. For this system, 
the variance of the fluctuations after appropriately normalizing, does 
not converge to the classical variance of the observable. However the 
classical variance can be recovered, after inserting some arithmetic cor- 
rection. Furthermore the limit distribution for the Hecke eigenforms 
is again not Gaussian (pending on conjectures of Keating and Snaith 



In this note we look at the matrix elements for quantized cat maps on 
T 4 , and the variance of their fluctuations. Here we find that for certain 
maps and observables, the variance has a different rate of decay than 
that predicted by the Feingold-Peres formula. 

Let A e Sp(4, Z) be a matrix with 4 distinct eigenvalues. For the 
dynamics to be ergodic we assume that A has no roots of unity for 
eigenvalues. We further assume that the vector space Q 4 decomposes 
into two rational orthogonal symplectic subspaces, invariant under the 
action of A. Denote by Z\,Z 2 , the lattices obtained by intersecting 
these subspaces with Z 4 . 

Analogously to the treatment of the case d — 1 we introduce Hecke 
operators, a family of commuting operators that commute with Un(A), 

and consider a basis of joint eigenstates {0j} J=1 , referred to as the 
Hecke basis. For any smooth observable / € C°°(T 4 ), denote the quan- 
tum variance in the Hecke basis by 



Theorem 1. Let f(x) = e 2mn ' x , n/fl. Then, as N — > oo through 
primes, the variance of these observables is given by 




CU). 



N 2 




+ 0(£) ifneZxUZ 2 , 
+ 0(-jkg) otherwise . 



This result is in contrast to the expected rate of predicted to 
hold for all observables in generic hyperbolic systems. 
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Remark 1.2. This result can be extended to any smooth observable 
/ e C°°(T 4 ), to state that S 2 (f, N) = ^ + O(^), if the Fourier co- 
efficients of / are supported outside the lattices ZiUZ 2 , and S 2 (f, N) = 
?P~ + O(jp) otherwise. The terms Vi(/), V 2 (f) can be expressed as 
quadratic forms in the Fourier coefficients of /, see ^2] for more details. 

After knowing the variance of the fluctuation, we can renormalize 
the matrix elements and consider the limit distribution. We present 
some numerical calculations of this distribution, for the two types of 
elementary observables (figures fip2l respectively). The numerical evi- 
dence suggests that when n G Z\ U Z 2 , the distribution of the normal- 
ized matrix elements is given by the semicircle law, where otherwise 
after appropriately normalizing, it corresponds to the product of two 
semicircle random variables. This result, is in agreement with the con- 
jecture presented in [H>] for the limit distribution for the cat map on 
T 2 . 

Furthermore, assuming that all matrix elements are roughly of the 
same order, the square root of the variance obtained here, gives an 
estimate of the rate of their decay. In figure E3 we present numeri- 
cal calculation for this decay, strongly supporting the validity of this 
estimate. 
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2. Background 

The full details for the cat map and it's quantization can be found 
in for one dimensional system, and in QjJ ^] for higher dimensions. 
We briefly review the setup for the quantization of the cat map on T 4 . 
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2.1. Classical dynamics. The classical dynamics are given by the 
iteration of a symplectic linear map A G Sp(4, Z). 

x = (P ^ G T 4 i-H. Ax. 

Given an observable / G C°°(T 4 ), the classical evolution defined by A 
is / l— / ° A. If A has no eigenvalues that are roots of unity then the 
classical dynamics is ergodic and mixing 



2.2. Quantum kinematics. For doing quantum mechanics on the 
torus, one takes Planck's constant to be 1/N, as the Hilbert space 
of states one takes Ti^ = L 2 ((Z/iVZ) 2 ), where the inner product is 
given by: 

y mod N 

For n = (ni,n2) G Z 4 (where Hi G Z 2 ), define elementary operators 
Tjy(n) acting on ip G Hn via: 

(2.1) T N (n)ip(y) = e N ( Ul ™ 2 )e N (n 2 ■ y)^(y + H\), 

where e^(x) = e^ 21 

For any smooth classical observable / G C°°(T 4 ) with Fourier ex- 
pansion f(x) = J2ftez 4 f(n)e 2mn ' x , its quantization is given by 

Op N (f) = E f(n)T N (n). 

Hez 4 

The main properties of the elementary operators TJv(n) are summa- 
rized in the following proposition. 

Proposition 2.1. For the operators Tjy(n) defined above: 

(1) ljv(n) only depends on n mod 2N. 

(2) Tjv(n) are unitary operators. 

(3) T7ie composition of two elementary operators is 

T N (m)T N (n) = e 2 N(^(rn,H))T N (H + m), 

where u(m,n) = ihx-n^—m-i-nx, is the symplectic inner product. 

(4) 7/n = (mod iV), tfien T^n) = (-l)^ 1 ^)/. 

These properties are easily derived from equation 12.11 Moreover, 
these properties uniquely characterize these operators, in the sense that 
if Tj^(n) are operators acting on an iV 2 -dimensional Hilbert space ful- 
filling the above properties, then they are unitarily equivalent to Tjy(n). 
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2.3. Quantum dynamics: For A G Sp(4, Z) which satisfies certain 
parity conditions (namely A = I (mod 2)), we assign unitary oper- 
ators Un(A), acting on L 2 ((Z/iVZ) 2 ) having the following important 
properties: 

(1) "Exact Egorov": For all observables / G C°°(T 4 ) 

UN{AY 1 Op N (f)U N {A) = Op N (f o A). 

(2) The quantization depends on A modulo 2N. 

(3) The quantization preserves commutation relations: 

if AB = BA (mod 2N) then U N (A)U N (B) = U N {B)U N {A). 

Note that if we multiply the operators by arbitrary phases, then all 
the above properties still hold. In fact the converse also holds, that 
is if Un(A) is an operator satisfying these properties, then U^(A) = 
e ia U N (A). 

Remark 2.1. Eventually, we are interested in the eigenfunctions of these 
operators, which are not affected by the choice of phase. Thus, in order 
to simplify the discussion we do not set explicit phases for the quantiza- 
tion. However, we note that one can make a specific choice of phases for 
which this quantization is multiplicative, Un(AB) = Un(A)Un(B). In 
particular this implies property (3) above, indeed holds for any choice 
of phases. 

Since we assume that N is odd, and the quantization factors through 
the group 

{B G Sp(4,Z/2iVZ)|£ = I (mod 2)} ~ Sp(4,Z/iVZ), 

we can define a (projective) representation of Sp(4, Z/iVZ): 
For any element B G Sp(4,Z/iVZ), let J? G Sp(4,Z/2iVZ) be the 
unique element s.t. B = B (mod N), B = I (mod 2), and identify 
Un{B) := Un(B). For the rest of this note we will use this convention 
for the quantization of Sp(4, Z/iVZ). 

2.4. Hecke eigenfunctions. Let A G Sp(4, Z), A = / (mod 2) and let 
C(N,A) be a maximal commutative subgroup of Sp(4, Z/iVZ), which 
includes A modulo N. 

Remark 2.2. In general, this group is not necessarily unique. However 
if all the eigenvalues of A mod iV are distinct then C(N,A) is unique. 
For an explicit construction see jT2l . 

Since the quantization preserves commutation relations, the opera- 
tors Un(B), B G C(N, A) form a family of commuting operators, called 
Hecke operators. Functions (p eTLn that are simultaneous eigenfunc- 
tions of all the Hecke operators are called Hecke eigenfunctions, and 
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a basis consisting of Hecke eigenfunctions is called a Hecke basis. By 
definition A G C(N,A), consequently any Hecke eigenfunction is in 
particular also an eigenfunction of Un(A). 

3. Reduction to lower dimension 

Let A G Sp(4, Z), A — I (mod 2), with 4 distinct eigenvalues. Fur- 
ther assume that the vector space Q 4 , decomposes into two (rational) 
symplectic subspaces, invariant under the action of A. In this case, the 
Hecke group C(N, A) = C±(N, A) x C 2 (N, A) is a direct product, each 
term Ci(N, A), can be identified with a lower dimensional Hecke group 
of a corresponding matrix in SL(2, Z/JVZ). Correspondingly, the Hecke 
operators and eigenfunctions are a tensor product of the appropriate 
lower dimensional Hecke operators and eigenfunctions. 

3.1. Reduction of Hecke group. For each invariant subspace, take 
a symplectic basis ti-.fi £ Q 4 ) that is Lo{e^e^) = cu(fi, fj) = and 
u(ei, fj) = 5ij. For a sufficiently large prime N, take ei, fi G (Z/iVZ) 4 , 
through reduction of and fi modulo iV respectively. This induces a 
decomposition of (Z/iVZ) 4 = E 1 © E 2 into two orthogonal symplectic 
subspaces invariant under the action of A mod N . 

For n G Z 4 , denote by (ni,n 2 ) G (Z/NZ) 2 x (Z/iVZ) 2 the restric- 
tion of n mod iV, to i?i © £2 in the symplectic basis. Since this is a 
symplectic decomposition then for any n, m G Z 4 

(3.1) uj(n, m) = uj{fii, fhi) + cj(n 2 , n 2 ) (mod AT). 

This decomposition induces a map 

i N : SL(2,Z/A^Z) x SL(2,Z/A^Z) Sp(4,Z/A^Z). 

If we denote by A G SL(2, Z/A^Z) the restriction of A (mod N) to ^ 
in the symplectic basis, then i^{A\, A2) = A (mod AT). 

Lemma 3.1. Let Ci(N,A) C SL(2,Z/N'Z) be maximal commutative 
subgroups which include A i; i — 1,2 respectively. Then the map i N 
defined above, maps the group Ci(N, A) x 2 (N, A) isomorphically onto 
C(N,A). 

Proof The map i N : Ci(AT, A) x C 2 (N,A) ^ C(N,A), is clearly in- 
jective. Thus, it is sufficient to show that this map is onto. Let 
B G Sp(4, Z/ATZ), be a matrix that commutes with A (mod N). We 
can assume A^ is large enough so that A (mod N) has 4 distinct eigen- 
values in F N 2 . Consequently, the spaces Ei are also invariant under the 
action of B. Thus if we denote by B { G SL(2, Z/NZ) the restriction of 
B to Ei, in the symplectic basis, then in(Bi, B 2 ) — B. □ 
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3.2. Quantization of Hecke group. Let Tjp(-), U^\-) be the quan- 
tized elementary observables and propagators for T 2 . For n G Z 4 , 
identify fii,n 2 defined above with elements of (Z/2iVZ) 2 by requiring 

(3.2) ni = (711,713), n 2 = (n 2 ,n 4 ) (mod 2). 

We can also identify B { G Ci(N,A), as elements of SL(2,Z/2iVZ) by 
requiring them to be congruent to / modulo 2. 

Proposition 3.2. There is a unitary mapping 

U : L 2 (Z/NZf -»■ L 2 (Z/iVZ) g> L 2 (Z/iVZ), 

such that: 

(1) For any n G Z 4 , 

UT N {n)U- 1 = T^^) ® T^\n 2 ). 

(2) For any B G C(iV, A), 

UU N {B)U- X = U$\B X ) ® U^\B 2 ). 

Proof. It is easily verified from l3.lll3.2l that Tjy(n) = T^\fii)®Tjp {n 2 ) 
obey the same relation as in proposition 12.11 Thus, from uniqueness 
they are unitarily equivalent. 

As for the second part, recall ujp(Bi) and U^\b 2 ) both satisfy the 
Egorov identity, and we showed that UT N (n)U" 1 = T^\ni) ®T^(n 2 )- 
Consequently, if we define U N (B) = U^ 1 U^\B 1 ) ® ujp(B 2 )U, then 
Un(B) satisfies the Egorov identity as well: 

U N {B)- l T N {n)U N {B) = T N {nB). 

Thus, from the uniqueness of the quantization 

UXJ N (B)U- X = US\B X ) <g> Ug\B 2 ). 

□ 

Remark 3.1. In fact, one can use this procedure to define the phases for 
the quantization of the Hecke group. Consequently, if the phases for the 
quantized two dimensional maps are chosen so that the quantization 
is multiplicative (as done in |15j). then the quantization of the Hecke 
group would be multiplicative as well. 



QUANTUM VARIANCE ON THE 4-DIMENSIONAL TORUS 



3.3. Hecke eigenfunctions. Let {<fij} N =1 , {$}* =1 e L 2 (Z/NZ), be 

joint eigenfunctions of all the Hecke operators U$(Bi), Bi G Ci(N, A) 
i — 1,2 respectively. Correspondingly, 0j lj2 = <8> 0| 2 ), are 

eigenfunctions of all the operators {t/jv(-B)|-B £ and hence 

{<f>ji,j 2 }j 1 j 2= i is a Hecke basis of L 2 (Z/NZ) 2 . In this basis, the matrix 
elements of the elementary observables Tjv(n), are given by 

(3.3) (T^fO^,^) = (^(n!)^,^)^^)^,^). 

Remark 3.2. The joint eigenspaces of all the operators 11$ (Bi), Bi G 
Ci(N,A), are one dimensional (except for the eigenspace correspond- 
ing to the trivial character) J2|- Correspondingly, any Hecke eigen- 
function is of the form defined above, except for Hecke eigenfunc- 
tions corresponding to the trivial character, which are of the form 
(f) = U- 1 (a(f) 1 1 ® 4>\ + b(f)\ <g> 2 ), where 4>\ G L 2 (Z/NZ) correspond to 
the trivial character, and \a\ 2 + \b\ 2 = 1. 

4. Variance of matrix elements 

We now turn to prove theorem^ We will concentrate on the elemen- 
tary observables f(x) = e 2nn ' x with corresponding quantum operators 
Op A r(/) = T N (n), and calculate their variance in the Hecke basis. 

Following the construction of the Hecke eigenfunctions and matrix 
elements described in the previous section, it is sufficient to understand 
the distribution of the matrix elements for the cat map on T 2 . In the 
following proposition, we summarize some results regarding the matrix 
elements of elementary observables in the Hecke basis on the 2-torus. 

Proposition 4.1. Let A G SL(2,Z/NZ) with 2 distinct eigenvalues. 
Let C(N,A) C SL(2,Z/NZ) be the corresponding Hecke group, and 
{4>i}f =1 G L 2 (Z/NZ) the corresponding Hecke eigenfunctions. For n G 
(Z/2NZ) 2 that is not an eigenvector of A (mod N): 

(1) The second moment of the corresponding operator is 

1 l/T'Wf-'U am 2 -] 1 ^ = (mod N) 

Nf^ l{ N W^'^'/l -\i + 0(^) n^O(modiV) 

(2) For cf), <f> Hecke eigenfunctions corresponding to the trivial char- 
acter: 

i^W^-u ,A| / 1 ^ = (mod N) 
\(T N (n)0,0)| = | m _ 1/2) __, Q (modAr) 

For proof we refer to ^H] • We now give the proof of theorem El 
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Proof. Let {(pj 1 ,j 2 }j 1 J2=1 be the Hecke basis constructed in the previous 
section. From the formula given for the matrix elements in 13.31 we can 
rewrite the quantum variance as 

N N 

S,{f,N) = 1 \(T%\n X )<t>Ul)\ 2 ^ E \^\n 2 W n Al)\ 2 . 

.71=1 J2=l 

Now if n G -Zi, then it is a linear combination of ex, fx, and hence 
for A" sufficiently large / 0, n 2 = mod iV. In the same way 
if n G Z 2 then fix = 0, n 2 / mod N, and if n ^ Zi U Z 2 then 
nx,n 2 7^ mod AT. Next, since A has no rational eigenvectors we can 
assume A" is sufficiently large so that fix,n 2 are not eigenvectors of 
Ax,A 2 respectively. The proof now follows directly from the first part 
of proposition 14.11 □ 

Remark 4.1. Note that if we take a different Hecke basis from the one 
constructed in the previous section, then we only change the elements 
in the sum corresponding to the trivial character. These elements, 
from the second part of proposition 14. 11 contribute O(-jb) (respectively 
0(jjz) if n G" Z\ U Z 2 ). Therefor this result holds for any Hecke basis. 

5. Discussion 

5.1. Product behavior. The main ingredient in the proof of theorem 
[U is the fact that the quantized system can be identified as a tensor 
product of two quantized two dimensional systems, and that there are 
observables whose quantization is trivial on one of the factors. 

This behavior is clearly to be expected if the classical four dimen- 
sional system is a product of two different two dimensional systems 
(e.g., considering the four dimensional torus product T 4 = T 2 x T 2 , 
and constructing an element A 6 Sp(4, Z) by taking two elements 
Ax,A 2 G Sp(2, Z) acting on each factor). In such a system, when look- 
ing at observables that are constant on one factor, the fluctuations of 
their quantization only come from the nonconstant factor, and hence 
the anomalous rate of decay. 

However, it is interesting to note that this behavior can also occur 
for systems that do not factor in terms of their classical dynamics. 
The condition for the classical dynamics to factor, is that there is a 
symplectic rational invariant decomposition such that the lattice Z\ © 
Z 2 = Z 4 (where Zx,Z 2 are as in theorem P). In some cases, even 
though the space Q 4 decomposes into two invariant rational symplectic 
subspaces, the lattice Z\ © Z 2 C Z 4 is a proper sublattice of finite index 
and the classical dynamics does not factor. Nevertheless, for any prime 
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iV the vector space (Z/iVZ) 4 decomposes in to two invariant symplectic 
subspaces, and the quantized system would still behave like that of a 
product of two systems. 

Remark 5.1. It is worth mentioning that even if the matrix A does not 
factor over the rationals, it does factor over a quadratic extension, and 
hence for half the primes it would factor modulo N. Consequently, for 
these primes we would still have product behavior, that is the quantized 
system is a tensor product of two quantized two dimensional systems. 
However, in this case there can be no classical observables whose quan- 
tization is trivial on one of the factors (for infinitely many primes). 

5.2. Limit distribution. The numerical evidence presented in figures 
11121 suggests that for the degenerate case (where the rate of decay is 
jj) the limit distribution follows the semicircle law, while otherwise 
it is given by the product of two random variables with semicircle 
distribution. 

In fact, assuming the validity of the Kurlberg-Rudnick conjecture 
regarding the limit distribution of matrix elements on the two dimen- 
sional torus jT5|, these results can be derived directly from the product 
behavior of the quantized system ()3.3j) . 

5.3. Rate of decay. The numerical evidence presented in figure H3 
imply a bound on the individual matrix elements: 



Given the Kurlberg-Rudnick rate conjecture (which is now a theorem 
|H]) the above bound (|5.1|) . can be derived from (|3.3|) as well. 



(5.1) 
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Figure 1. Scaled distribution of normalized matrix el- 
ements Fi{n) = y/N(T N (n)<j>i,<l>i) , n = (1,0,0,0) e Z x 
for N = 1613. Plotted together with a semi circle distri- 
bution y(t) = - t 2 . 
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-4-3-2-1 1 2 3 



Figure 2. Scaled distribution of normalized matrix ele- 
ments Fi{n) = N{T N (n)<Pi, , n= (1,1,0,0) £ Z X \JZ 2 
for N = 1613. Plotted together the corresponding dis- 
tribution y(t) = £ f* ^(4 - s 2)(4 - (f ) 2 )f 
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Figure 3. The values of maxj |(T/v(ra)</>i, (f>i)\ for n = 
(1,0, 0,0), (1,1, 0,0) and values of N between 100 and 
1000. Plotted together with the predicted decay of 
2iV~ a5 , 4iV~ 1 respectively, on a log-log scale. 
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